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TITLE:  STUDIES  ON  THE  COMBUSTION  EFFICIENCY  OF 

LIQUID  FUEL  ROCKETS  WITH  VARIABLE  THRUST 

AUTHOR:  Li  Xiaobin,  Chen  Qizhi,  Chen  Yunqin 

SUMMARY 

This  article  carries  out  studies  on  the  stability  characteristics 
of  liquid  fuel  rocket  engines  within  a  relatively  large  range  of 
variations.  The  research  work  was  carried  out  centering  around  the 
speed  and  efficiency  characteristics  of  combustion  chambers.  This 
article  sets  out  the  results  of  relevent  experiments  and  theoretical 
calculations.  From  the  angle  of  the  combustion  processes  within 
combustion  chambers/  there  were  clear  indications  of  the  status  of  the 
effects  of  various  processes  (for  example/  evaporation  or 
vaporization,  mixing,  chemical  reactions,  etc.)  following  along  with 
changes  in  the  operatinal  status  of  variable  thrust  engine  combustion 
efficiency . 

SYMBOL  EXPLANATIONS 

2 

A  combustion  chamber  cross  section  area  cm  ;  m  propellant 
c  2  ^ 

mass  flow  kg/s;  A.  jet  tube  throat  area  cm  ;  ni  ^  liquid 

t  vap 

drop  evaporation  rate  kg/s;  C  constant  pressure  specific  heat 
o  P  3  * 

J/kg *  K;  N  liquid  drop  number  in  unit  volume  1/cm  ;  C 

2 

characteristic  speed  m/s;  p  pressure  N/m  ;  D  diffussion 

2  ... 
coefficient  m  /s ;  Q  heat  added  to  liquid  drop  m  unit  time  by 

gases  J/s;  E  decomposition  reaction  activation  energy  J/kgmol ;  P 

effective  or  applied  force  of  gases  on  liquid  drop  kg;  Q 

c  2 

decomposition  reaction  heat  J/kg;  g  strength  of  flow  kg/s  m  ;  $ 

2 

heat  flow  strength  vector  J/s*m  ;  gn  initial  propellant 

O 

evapor i zation  or  vaporization  fraction  ;  gQ  heat  liquid  drop  heat 
exchange  rate  J/s;  H  total  liquid  drop  unit  mass  enthalpy  J/kg;  R 
universal  gas  constant  J/kg  o°K;  h  degree  of  opening  of  jet 
injector  needle  valves  mm;  r  liquid  drop  radius  m;  total 

gas  unit  mass  enthalpy  or  total  heat  ;  Re  Reynolds  number;  H 

§ 

total  enthlpy  of  liquid  drop  unit  mass  J/kg;  T  temperature  K; 

K  component  ratio;  fi  gas  speed  m/s;  Kq  decomposition  reaction 
constant  1/s;  P  liquid  drop  speed  m/s;  L  characteristic 
length  m;  a  surplus  oxygen  coefficient;  m  liquid  drop  mass  kg;  y 
specific  heat  ratio;  M  molecular  weight  kg/gkmol;  6  liquid  drop 
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distribution  standard  deviation;  n  efficiency;  p-  liquid  drop  mass  in 
unit  volume  kg/m3;  A  heat  conduction  coefficient  J/s  »m  *°K  ;*>*  gas 

density  kg/m3;  it  viscosity  k-s/m2;  x  (illegible)  lii<3uid' 
drop  density  kg/m3;  v  decomposition  reaction  series;  ?  friction 

force  tensor  kg/m3. 

SUBSCRIPT  EXPLANATION 

a  adiabatic;  min  minimum;  c  combustion  chamber  environment;  mix 
mixture;  f  fuel;  n  nth  dimension  set;  F  decomposition  flame 
peak;  o  oxidizer-initial  value,  basic  value;  g  gas;  s  liquid 
drop  surface  area;  i  ith  root  flow  tube;  t  total  volume;  ing  jet 
injection;  th  theoretical  value;  j  jth  type  component;  vap 
vaporization;  1  liquid  drop;  exp  experimental;  max  maximum 
value;  cool  cool  or  cold  test 

I.  INTRODUCTION 

Following  along  with  the  daily  development  of  space  flight 
technology,  the  requirement  for  variable  thrust  liquid  fuel  rocket 
engines  is  becoming  more  and  more  pressing.  In  spacecraft  attitude 
control  for  such  areas  as  orbital  control  and  planetary  and  stellar 
surface  probes,  variable  thrust  liquid  fuel  rocket  engine 
applications,  in  all  cases,  are  advantageous  or  indispensible .  As  far 
as  the  wide  range  of  changes  in  the  operating  conditions  of  variable 
thrust  liquid  fuel  rocket  engines  is  concerned,  in  terms  of  energy 
transfer  processes  within  combustion  chambers,  they  present  even 
higher  requirements.  Because  of  this,  how  to  guarantee,  within  the 
whole  range  of  thrust  changes,  the  achieving  of  combustion 
efficiencies  as  high  as  possible  and  as  stable  as  possible  is  a 
research  problem  wihich  possesses  practical  weapons.  In  experimental 
research,  one  discovers  that  changes  in  the  operational  configuration 
of  variable  thrust  liquid  fuel  rocket  engines  have  very  great  effects 
on  the  combustion  efficiencies  produced.  A  5:1  thrust  variation  ratio 
liquid  fuel  rocket  engine,  within  the  entire  operational  flow 
adjustment  range,  has  changes  in  combustion  efficiency  which  are 
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capable  of  dropping  from  90%  of  design  configuration  to  a  low  of 
approximately  60%.  We  used  a  125kg  level  experimental  engine  and 
carried  out  theoretical  and  experimental  research  on  the  combustion 
efficiencies  in  variable  thrust  operations. 

II.  THEORETICAL  MODELS 

As  far  as  the  setting  up  of  numerical  models  of  combustion 
..processes  within  combustion  chambers  is  concerned,  using  analytic 
methods  in  order  to  describe  the  interior  energy  transfer  processes  of 
combustion  chambers  is  a  problem  which  people  have  been  probing  right 
along.  Since  1946,  scores  of  models  have  been  presented.  As  far  as 
the  utilization  of  numerical  models  is  concerned,  it  is  possible  to 
understand  in  relative  depth  the  degree  of  the  effects  of  changes  in 
the  various  operational  parameters,  within  combustion  chambers,  when 
there  are  thrust  variations,  as  these  changes  effect  all  the  various 
energy  transfer  processes.  As  a  result  of  this,  one  understands 
control  factors  and  the  change  configurations  for  energy  transfer 
during  thrust  changes  and  deeply  comprehends  the  causes  for  changes  in 
combustion  efficiency  which  follow  changes  in  thrust. 

We  know  that  energy  transfer  processes  within  combustion  chambers 
are  composed  of  a  number  of  complicated  subprocesses.  Speaking  in 
terms  of  simple  generalization,  they  mainly  accompany  the  exchange 
processes  associated  with  high  speed  chemical  reactions  in  two-phase 
gas-liquid  fields.  Jet  injection  vaporization  or  atomization 
processes  are  relatively  complicated.  It  is  difficult  to  use  analytic 
methods  to  describe  them.  At  the  present  time,  one  is  still  dependent 
on  experimentation.  The  speed  of  gas  phase  chemical  reaction 
processes  is  relatively  fast.  It  is  not  easy  to  turn  them  into  energy 
trasfer  control  processes.  Because  of  this,  models  take  as  a 
condition  the  complete  disintegration  of  liquid  drops  in  vaporization 
or  atomization.  In  conjunction  with  this,  it  is  recognized  that  gas 
phase  chemical  reactions  are  instantaneous.  These  models  principally 
describe  flow  fields  aft  of  jet  injection  vaporization  or  atomization 
zones.  In  conjunction  with  this,  they  lay  emphasis  on  liquid  drop 
evaporation  or  vaporization  diffusion  processes. 
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If  one  does  not  consider  liquid  drop  mutual  collisions  and 
secondary  rupture,  as  far  as  liquid  drop  is  concerned,  one  has  the 
equations  set  out  below  for  liquid  drop  quantitative  magnitude,  mass, 
momentum,  and  energy 


(1) 


d"dT~!~  *  V ' iP '  9  ‘  V '  3  "  A'  ‘ * :  ~ 'V ;  ,il' • ' r '  (3) 


In  these 


p;-.v>; 

As  far  as  gases  are  concerned,  it  is  also  possible  to  set  out  the 
equations  below  for  mass,  momentum,  energy,  as  well  as  state 
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If  one  assumes  that  flow  fields  are  one  dimensional,  and,  at  the 
same  time,  only  considers  stable  state  processes,  not  considering  the 
non-homogeneous  natures  of  heat  exchange  friction,  and  combustion 
chamber  pressures,  one  uses  localized  propellant  components  in 
accordance  with  thermodynamic  calculations  to  determine  combustion  gas 
temperature.  In  this  way,  equations  (1)  to  (8)  described  above  are 
capable  of  being  simplified  into  the  equations  set  out  below 
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In  these,  K  is  the  local  component  ratio 
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uv  is  the  terminal  velocity  that  hot  gases  are  capable  of 
reaching 
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In  equation  set  (9),  for  the  liquid  drop  gas  phase  applied  force 
relationship  quantity  F,  one  opts  for  the  expression  below 
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In  this,  as  far  as  the  liquid  drop  evaporation  or  vaporization 
amount  m  is  concerned,  considering  the  heat  decomposition 

V  3  p 

reaction,  it  is  determined  by  the  several  equilibrium  relationships 
specified  below 
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As  far  as  the  heat  exchange  quantity  0  is  concerned,  H  is 

O 

related  as  shown  below 
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Opting  for  the  use  of  reasonable  boundary  conditions  and  special 
thermal  characteristic  calculation  methods  makes  it  possible  to  solve 
equation  set  (9)  for  numerical  values.  As  a  result  of  this,  it  is 
possible  to  solve  for  and  obtain  propellant  evaporation  or 
vaporization  processes  along  the  direction  of  the  axis  of  the 
combustion  chamber.  In  conjunction  with  that,  finally,  one  obtains 
evaporation  or  vaporization  rates  within  combustion  chambers. 

III.  EXPERIMENTAL  EQUIPMENT  AND  TEST  RESULTS 
Fig.l  shows  a  schematic  diagram  of  the  system  used  for 
experimental  research.  The  system  as  a  whole  is  composed  of  4 
subsystems.  From  left  to  right,  respectively,  are  the  oxidation  agent 
transport  system,  the  control  water  supply  system,  the  cooling  water 
control  system,  and  the  fuel  transport  system.  Experimental  engines 
opt  for  the  use  of  coaxial  splash  or  spray  ring  type  jet  injector 
devices.  Thrust  adjustment  is  effected  through  the  degree  of  opening 
of  jet  injection  paths  by  which  a  moveable  needle  valve  simultaneously 
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Table  1  Fig.l  Symbol  Explanations  (1)  Gas  Source  (2)  Pressure 
Reduction  Device  (3)  Oxidizer  Storage  Tank  (4)  Gas  Release  Valve 
(5)  Injection/Exhaust  Valve  (6)  Main  Oxidizer  Valve  (7)  Filter 
Net  (8)  Flow  Meter  (9)  Pneumatic  Valve  (10)  Blow  Off  Valve  (11) 
Fuel  Storage  Tank  (12)  Main  Fuel  Valve  (13)  Control  Water  Storage 
Tank  (14)  Main  Control  Water  Valve  (15)  Electromagnetic  Entry 
Valve  (16)  Electromagnetic  Exit  Valve  (17)  Cooling  Water  Storage 
Tank  (18)  Main  Cooling  Water  Valve  (19)  Jet  Injector  Device  (20) 
Thrust  Chamber 


3) 

<3>£i£ 


12S~*2$kf 
lfl~2  kg/cm't*) 
0.5~].0tn 
70  m  m 


32mm 


Table  2  Main  Thrust  Chamber  Design  Parameters  (1)  propellant: 
nitric  acid  27/nonsymetr ical  dime|hylhydrazine  (2)  thrust:  125-25kg 
(3)  chamber  pressure:  10-2kg/ciri  (absolute)  (4)  characteristic 
length:  0.5-1.0m  (5)  combustion  chamber  diameter.  70mm  (6)  jet 

tube  throat  diameter:  32mm 

The  experimental  engine  carries  out  large  amounts  of  heat 
testing.  Fig. 2  and  Fig. 3  show  partial  results.  In  the  Fig.'s 
corresponding  pressure  is  ~  p  /p  .  Fig. 2  opts  for  the  use 

of  a  water  cooled  thrust  chamber,  and  one  time  test  results  for 

* 

variable  thrust  engines  with  characteristic  length  L  =  0.5mm. 

These  clearly  show  engine  combustion  efficiency  v c  *,  a  needle  valve 

opening  degree  of  h,  and  a  jet  injection  pressure  drop  of  .a.  p.  , , 

l ,  t 

following  along  with  changes  in  operational  status.  Fig. 3  is  engine 
variable  thrust  operation  results  associated  with  opting  for  the  use 

it 

of  water  cooled  thrust  chambers,  L  =1.0m. 

From  operational  results,  it  is  possible  to  see  that: 

1)  in  the  process  of  adjusting  engine  thrust  from  25kg-125kg, 
combustion  chamber  pressure  p  ,  jet  injection  pressure  drop  A.p., 

and  the  needle  valve  degree  of  opening  h  all  go  through  broad  changes, 

2)  combustion  chamber  characteristic  speed  efficiencies,  when 

adjusted  within  a  range  of  25-125kg,  show  broad  changes. 

* 

This  is  particularly  obvious  in  the  case  of  L  =0 . 5m  combustion 
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chambers.  When  >c>0.8,  efficiency  *  is  capable  of  reaching 
over  90%,  However,  when  >,-<0.5,  7  *  drops  below  80%.  At  times 
of  minimum  thrust,  7  changes  are  the  most  abrupt. 

3)  As  far  as  the  rules  or  patterns  described  above  for  changes 

in  * _ p  are  concerned,  it  goes  without  saying  that  they  are 

c 

consistent  whether  for  water  cooled  thrust  chambers  or  combustion 
corrosion  thrust  chambers. 

4)  When  there  are  changes  in  combustion  chamber  characteristic 

* 

lengths  L  ,  there  is  a  relatively  large  effect  on  • tj 

c  c 

relationship  curves.  In  small  operational  configurations,  this 

* 

influence  is  most  obvious.  For  example,  when  L  increases  from  0.5m 

to  1.0m,  7,  *  values  corresponding  to  small  p  values  show 

c  c 

relatively  large  increases. 

IV.  COMBUSTION  EFFICIENCY  ANALYSIS 

In  the  general  run  of  liquid  fuel  rocket  engines,  due  to  the 
majority  of  jet  injection  device's  being  composed  of  multiple 
individual  jet  apertures,  aft  of  the  jet  injection  area,  the  radial 
direction  mixing  is  relatively  weak.  Because  of  this,  it  is  possible 
to  take  the  combustion  and  flows  within  combustion  chambers  and  divide 
them  into  a  number  of  flow  tubes  parallel  to  the  axis  line  for 
purposes  of  consideration.  Each  flow  tube  operates  like  a  little 
independent  combustinon  chamber.  Because  of  this  the  combustion 
efficiency  7  c*  is  capable  of  being  expressed  as  the  sum  of  the 
effects  of  the  various  flow  tubes: 


7c*  = 


• ,  "" 


1 f  one  carries  out  transformations  of  the  equation  above,  connecting 
it  to  the  initial  jet  input  component  ratio  distribution,  one  then  has 

« 

=  7c*.*r«#,7c*.i»i. 
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From  this,  we  will  take  combustion  efficiency  and  divide  it  into 

two  parts.  One  is  n. c*  vap»  it  is  related  to  propellant 

evaporation  or  vaporization.  It  is  called  vaporization  efficiency. 

The  second  is  rj,  *  .  .  It  is  related  to  propellant  mixing.  Its 

c,nux  3 

name  is  mixing  or  mixture  efficiency. 

1.  Mixing  Efficiency 

From  analysis  one  learns  that,  in  the  process  of  engine  thrust 
variations,  the  phenomenon  of  propellant  cavitation  in  tubes  will  not 
-occur.  The  "separation"  phenomena  of  jet  injector  ram  atomization  or 
vaporization  (that  is,  RSS  phenomena)  will  also  not  be  severely 
present.  Because  of  this,  in  our  analysis,  we  used  non-reactive 
simulation  media  cold  flow  tests  in  order  to  precisely  determine  rjf 
c*mix’  It:  to  a  definite  degree,  reflect  actual  propellant 

jet  mist  mixture  configurations.  Considering  ram  atomization 
processes  to  be  controlled  by  the  geometrical  dimensions  and  momentum 
of  jets,  in  tests,  one  opts  for  the  use  of  similar  conditions  in  which 
jet  injector  geometrical  dimensions  are  equal  and  momentum  is  equal. 

As  far  as  flow  strength  distributions  of  jet  vaporization  fields  and 
mixture  ratios  are  concerned,  precise  measurements  were  carried  out, 
and,  in  conjunction  with  that,  mixture  efficiencies  were  calculated 
out.  Fig. 4  and  Fig. 5  show  partial  experimental  results. 


Fig. 4  Flow  Strength  and  Mixture  Ratio  Distributions  Under 
Non-Specif ied  Conditions  (1)  Direction 
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Fig. 4  shows  jet  injector  flow  field  and  mixture  ratio 
distributions  for  small  thrusts.  At  this  time,  because  the  membrane 
thickness  of  jet  injection  jets  is  only  0.058mm,  as  a  consequence,  the 
jet  injection  speed  is  also  very  large.  Besides  this,  due  to  the  fact 
that,  with  small  thrusts,  the  degree  of  opening  of  jet  injector 
devices  is  very  small,  the  jet  display  for  the  smoothness  and 
geometical  position  tolerances  of  parts  is  very  sensitive.  These 
factors  make  jet  flow  membranes  or  films  very  unstable.  After  the 
jets  come  in,  they  then  split  up  to  become  thread  shaped  bundles  of 
liquid  even  causing  jets  to  show  the  appearance  of  clear  deflection. 

As  a  result  of  this,  the  effects  of  ram  mixing  between  two  types  of 
components  are  weakened  and  reduced  causing  flow  strength  and  mixing 
ratio  distributions  to  be  quite  uneven  .  In  Fig. 4,  the  displayed  flow 
strength  peak  values  and  peak  values  for  mixture  ratio  distributions 
are  precisely  reflections  of  the  effects  discussed  above.  From  this; 
it  is  clearly  demonstrated  that  the  results  of  propellant  mixing,  in 
small  operational  configurations,  are  bad. 

What  Fig. 5  shows  is  the  status  for  engine  spray  flow  strength  and 
mixing  ratio  distributions  with  specified  or  rated  thrusts.  At  this 
time,  jet  injector  jet  fluid  flims  or  liquid  membranes  are  relatively 
thick.  Jet  injection  speed  is  relatively  small.  The  influence  of  the 
processing  quality  of  jet  injection  devices  is  correspondingly  weak. 

As  a  result,  jet  liquid  spray  films  or  membranes  are  stable.  Results 
for  jet  injection  component  fluid  films  or  liquid  membranes  going 
through  ram  and  splash  or  spray  mixing  are  relatively  adequate.  The 
broken  line  obtained  from  experiments  for  high  operating 
configurations  are  also  relatively  round  and  smooth. 

On  the  basis  of  the  definitions  on  the  front  page,  and,  making 
use  of  cold  flow  test  results,  calculations  were  carried  out  for 
mixing  efficiencies  with  variable  operational  configurations.  As  to 
the  fact  that  sampling  was  not  carried  out  on  flow  fields  as  a  whole, 
a  choice  was  made  to  use  the  formula  below  in  order  to  calculate 
mixing  efficiencies: 


I 

2  £,C*(a,„) 


i3 


68 


(16) 


Fig. 5  Flow  Strength  and  Mixing  Ratio  Distributions  Under  Specified  or 
Rated  Conditions  (1)  Direction 


As  far  as  consideration  of  the  structural  symmetry  of  ring  shaped 
jet  injector  devices  which  were  studied  is  concerned,  in  tests, 
sampling  was  only  done  on  the  radial  directions  of  the  four  quadrants 
of  jet  spray  fields.  In  order  to  guarantee  precision,  eacn  set  of 
experiments  was  replicated. 


ft(m  m) 

1  0.255 

!  0.241 

!  0.234 

.'  0.J51  ! 

0.09 

0,(kg/cm:.) 

10.0 

7.1 

!  5-* 

3.8 

1.9 

*.(kg/J) 

0.399 

0.276 

I  0.2SI 

0.1*7 

0.090 

*  c  (kg/*) 

0.155 

;  0.147 

0.143 

0.016 

0.050 

*7e  *»«■!  Ik 

1.007 

j  1.001 

1.001  j 

0.959  j 

0.955 

Table  3  Mixture  Efficiency  Calculation  Results 
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Cold  test  parameters  are  designed  for  conditions  which  simulate  a 
certain  hot  test  condition.  Test  conditions  and  corresponding  cold 

test  results  are  set  out  in  Table  3.  In  the  table,  the  appearance  of 

* 

mixture  efficiencies  >  1  is  due  to  the  fact  that  C  (a.)  values  are 
*  ^ 
greater  than  C  (a^nj)  values,  which  partially  correspond  to  the 

surplus  oxygen  coefficient  a^  >  a^nj*  These  compensate  for  the 

partial  losses  of  a^  <  a^nj*  At  the  same  time,  it  is  also  due  to 

the  existance  of  differences  in  measurements.  From  Table  3,  it  is 

possible  to  see  that,  in  a  change  in  operational  status  or 

-configuration  from  high  to  low,  the  thinner  jet  injection  liquid 

membranes  or  fluid  films  (h)  are,  the  more  various  types  of  factors 

influencing  mixing  follow  this  and  aggravate  things.  This  leads  to  a 

drop  in  jet  device  mixture  efficiency.  Data  clearly  shows  that  when 

f.  >  0.6,  mixing  is  not  perfect,  but  losses  are  relatively  small. 

When  fc  <  0.6,  mixing  is  not  perfect,  and  losses  gradually  increase. 

The  lower  the  operational  configuration  is,  the  smaller  the  mixing 

efficiency',;^*,  mix,th  is. 

In  this,  mixing  efficiency  is  completely  dependent  on  the  results 
of  cold  tests  for  its  calculation.  In  actual  engine  operation,  due  to 
the  existence  of  combustion  reactions,  actual  mixing  states  and  cold 
flow  tests  will  show  some  differences.  Moreover,  combustion  chamber 
processes,  on  the  basis  of  flow  tube  division,  will  also  have  definite 
differences  from  the  actual  situation.  Because  of  this,  the  mixing 
efficiencies  obtained  by  calculation  will  not  completely  match  up  with 
the  actual  situation.  However,  cold  flow  tests  clearly  reflect  the 
important  fact  that  propellant  mixtures  in  combustion  chambers  follow 
changes  in  operational  configuration.  As  a  result,  they  clearly  show 
that  the  bad  turn  in  mixing  processes  is  an  important  factor  in  the 
lowering  of  combustion  efficiencies  under  low  operational 
conf igurations . 

2.  Vaporization  Efficiency 

As  far  as  evaporation  or  vaporization  efficiencies  are  concerned, 
we  opted  for  the  use  of  combustion  models  which  already  have 
vaporization  in  order  to  do  calculations.  It  is  already  known  that, 
in  the  process  of  variable  thrust  operations,  as  far  as  changes  in  the 
operating  conditions  for  engine  jet  injectors  and  combustion  chambers 

are  concerned,  one  sees  the  occurence  of  wide  changes.  Following 
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operational  configuration  from  high  to  low,  jet  injection  speed 
increases  continually.  Moreover,  the  degree  of  opening  of  jet 
injector  devices  decreases  continually.  These  changes  necessarily 
produce  effects  on  the  status  of  vaporized  or  atomized  propellent 
sprayed  into  the  interior  of  the  combustion  chambers. 

Detailed  analysis  of  vaporization  or  misting  processes  in 
experimental  engines  makes  it  possible  to  know  that  mutual  collision 
effects  associated  with  dual  component  fluid  films  or  liquid  membranes 
occupies  a  leading  position.  In  reference  [4],  the  liquid  drops 
produced  at  two  mutual  collision  type  jet  nozzles  are  given  by 


Dm  — 


i/;-“ 


From  the  similarity  of  atomization  processes  and  control  mechanisms, 
it  is  possible  to  recognize  in  an  approximate  way  that  there  is,  in 
our  experimental  engines,  the  relationship  below 


r  _  r  ,  -  A  g / 

'  T  ■!•*./  ~  **  y  «  1 7  %  j  , 


(17) 


In  calculations,  one  opts  for  the  use  of  equation  (17).  In  the 
equation,  b^  and  bQ,  respectively,  are  the  thicknesses  of 
propellant  liquid  membranes  or  films  and  oxidizer  liquid  membranes  or 
films.  A,  then,  is  indirectly  determined  utilizing  experimental 
results  for  a  certain  operational  configuration.  With  regard  to 
liquid  mist  distribution  parameters,  on  the  basis  of  reference  [3], 
for  the  liquid  drop  standard  deviation  <5  and  spectrum  width  Z,  one 
picks  <5  =  0.82  and  Z  =  3.0. 

In  evapor i zation  or  vaporization  efficiency  calculation  models, 

initial  spray  or  jet  atomization  zones  are  described  with  simplified 

models.  The  key  parameters  for  these  areas  include  jet  injection  jet 

disintegration  length,  L  ,  and  propellant  evaporation  or 

sp 

vaporization  amount  fraction  inegible )  *  From  t*16  3et  mist  or 
vapor  status  in  experimental  engines  discussed  previously  and  its 
special  features,  it  is  possible  to  see  that,  with  regard  to  the 
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non-uniform  characteristics  which  exist  in  conical  jet  atomization 
fields  as  well  as  they  themselves,  they  cause  the  radial  direction 
exchanges  associated  with  jet  injection  misting  zones  in  the 
combustion  processes  to  be  relatively  strong  and  violent.  As  a 
result,  reflux  and  turbulent  flow  strengths  are  all  relatively  large. 
This  will  strengthen  the  mixing  of  propellants,  and,  at  the  same  time, 
shorten  the  jet  decomposition  length.  When  operational  configurations 
change,  changes  in  combustion  chamber  conditions  will  have  relatively 
.great  effects  on  jet  decomposition  lengths.  When  thrusts  go  from 
large  to  small,  jet  speeds  go  from  small  to  large.  Under  small 
operational  configurations,  liquid  membrane  or  fluid  film 
disintegration  into  liquid  threads  will  cause  mutual  colllision  mixing 
results  between  components  to  vary.  At  the  same  time,  because  chamber 
pressures  drop,  this  causes  combustion  gas  densities  to  be  reduced. 
Liquid  phase  reactions  are  attenuated,  jet  injection  atomization  or 
misting  zone  turbulent  flow  strengths  are  reduced,  and  other  similar 
factors.  This  will  make  jet  decomposition  lengths  correspondingly 
longer . 

In  another  area,  when  one  uses  natural  propellants,  the  liquid 
phase  reactions  must  be  given  consideration.  In  calculation  models, 
this  partial  liquid  phase  reaction  is  expressed  by  gD.  From  the 

D 

mechanism  of  liquid  phase  reactions,  it  is  possible  to  know  that 
liquid  phase  reaction  speeds  and  combustion  chamber  pressures  are 
related  in  a  direct  proportion.  At  the  same  time,  the  level  of  liquid 
phase  reactions  is  also  in  direct  proportion  with  the  size  of  the 
areas  of  the  liquid  phase  contact  between  propellant  components. 
Following  along  with  the  alterations  in  operational  configuration 
described  above,  all  factors  are  in  flux.  Because  of  this,  the 
determined  g^  will  also  change  along  with  the  operational 
configuration.  In  high  operational  configurations,  high  combustion 
chamber  pressures  create  relatively  high  liquid  phase  reaction  speeds. 
At  the  same  time,  the  numerical  values  for  propellant  flow  strength 
are  large.  Moreover,  distributions  are  relatively  concentrated, 
mixing  is  relatively  perfect,  and  total  contact  area  between 
propellants  is  relatively  large.  Because  of  this,  in  high  operational 
configurations,  initial  evaporation  or  vaporization  amounts  for 
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propellants  are  relatively  large.  Following  drops  in  thrust,  initial 
evaporation  or  vaporization  amounts  are  reduced.  As  a  first 
approximation ,  take  L  an<3  g g*  Following  changes  in  operational 
configuration,  both  form  linear  relationships  with  chamber  pressures 
Pc ,  that  is, 

f./=4.0-0.248(/>e-i.J73) 

9*  =  0.05  +  0.0310(/C  —  1.873) 

In  these  equations,  the  various  constant  values  are  based  on  the 
general  range  of  gR  and  L  given  in  the  relevent  references.  At 
the  same  time,  consideration  is  also  given  to  the  actual  situation  of 
the  experimental  engine  in  question  as  well  as  consideration  being 
given  to  the  status  of  burn  corrosion  on  the  inner  walls  after  hot 
tests  of  combustion  corrosion  combustion  chambers  for  precise 
determinations . 

Going  through  the  analysis  and  processing  above,  it  is  then 
possible  to  carry  out  calculations  for  evaporation  or  vaporization 
efficiencies.  Calculation  conditions  and  mixture  efficiency  analysis 
conditions  are  consistent  with  each  other.  Actual  conditions  and 
results  are  set  out  in  Table  4. 

From  Table  4  it  is  possible  to  see  that,  following  along  with 
reductions  in  thrust,  the  vaporization  efficiency  of  experimental 
engines  drop. 

In  the  vaporization  efficiency  calculations  discussed  above,  they 
are  carried  out  on  the  basis  of  one  dimensional  flow  tube  theory.  No 
consideration  is  given  to  the  effects  of  non-uniform  spray  area 
propellant  distributions  on  evaporation  or  vaporization  efficiencies. 
This  influence,  under  low  operational  configurations,  will  not  be 
negligible.  Because  of  this,  the  evaporation  or  vaporization 
efficiencies  obtained  by  calculation  follow  changes  in  operational 
conf iguration  and  show  a  certain  disparity  with  the  actual  situation. 
However,  in  terms  of  mechanism,  they  illucidate  the  trends  in 
evaporation  or  vaporization  efficiencies  as  they  follow  changes  in 

operational  configuration.  They  clearly  show  that,  when  under  low 

* 

operational  configurations,  L  =0.5m  engines  have  relatively  large 
effects  due  to  imperfect  evaporation  o r  vapor ization  lowering 
combustion  efficiencies. 
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Table  4  Evaporation  or  Vaporization  Efficiency  Calculation  Results 


3.  Comprehensive  Analysis  of  Characteristic  Speed  Efficiency 


On  the  foundation  of  solutions  for  »)4c*,mix,th  and  v 
c* , vap , th  obtained  in  cold  flow  tests  and  theoretical  calculations , 
from  f?c*,th  =  "17  c*  ,mix  ,th  '(illegible)  irc*,vap,t,  it  is 
possible  to  obtain  the  characteristic  speed  efficiency  of  theoretical 
analysis  v  *,th.  In  Fig. 6,  one  takes  changes  in  the  theoretical 
values  17  *,th  and  the  test  values  '  rj  *,exp  and  makes  a 
comparison.  In  the  Fig.,  = 

From  Fig. 6  it  is  possible  to  see  that.r?c*,th  and  n  ^*,exp 
are  consistent  in  the  overall  rules  or  patterns  they  obey  in  following 
changes  in  operational  configuration.  Proceeding  to  the  next  step, 
this  explains  this  article's  being  directed  at  the  analytical  models 
discussed  above  which  were  set  up  for  variable  thrust  rocket  engines. 
In  terms  of  mechanisms,  they  reflect  factors  influencing  combustion 
processes  inside  actual  engines. 

In  Fig. 6,  one  sees  reflected  the  fact  that  i)i  c*  ,th  and  9 
*,exp  have  associated  with  them  patterns  or  rules  in  their  changes 

v 

which,  in  terms  of  actual  trends,  still  show  the  existence  of 
differences.  The  reasons  are:  1)  when  calculating  17  c*,vap,th, 
there  is  no  consideration  given  to  the  influence  of  non-uniform  mixing 
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distributions  in  spray  zones.  This  influence,  during  low  operational 


Fig. 6  Comparison  of  *  Values  from  Theoretical  Analysis  and 
Experimentation  (1)  Theoretical  Analysis  (2)  Experimentation 


Fig. 7  A  Comparison  of  Efficiencies  with  Different  Combustion  Chamber 
Lengths  (I)  Solid  Symbol  L  =1.0m  Hollow  Symbol  L  =0.5m  (1) 

Test 


configurations,  is  particularly  obvious.  2)  When  calculating  rj’ 
*vap,th,  assume  that  the  gas  phase  chemical  reactions  are  completed 

V 

instantaneously,  and  they  do  not  produce  effects  on  combustion 
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efficiencies.  Moreover,  when  actual  chamber  pressures  are  relatively 
low,  chemical  reaction  speeds  will  drop.  Actual  combustion  gas 
temperatures  will  be  lower  than  equilibrium  temperatures.  As  a  result 
of  this,  effects  will  be  produced  on  evaporation  or  vaporization 
speed.  3)  '  i),  c*mix,th  is  obtained  based  on  cold  flow 
experimentation.  It  shows  certain  discrepancies  with  actual  hot  test 
states.  4)  Under  various  operational  configurations,  actual 
measurements  of  liquid  jet  sprays  were  not  carried  out.  Generalized 
applications  were  made  of  the  test  results  of  previous  people.  It  is 
difficult  to  precisely  reflect  the  actual  rules  or  patterns  for  the 
jet  injector  devices  in  question.  5)  In  calculations,  we  ignored 
heat  transmission  and  friction  losses. 

In  order  to  carry  a  step  further  checks  on  the  level  of  influence 
on  energy  transfer  efficiencies  of  evaporation  or  vaporization  and 
mixing  when  testing  engines  under  varied  operational  configurations, 
one  opts  for  the  use  of  different  lengths  of  combustion  chamber  to 

carry  out  hot  tests.  Experimental  results  are  shown  in  Fig. 7.  When 

* 

L  increases  from  0.5m  to  1.0m,  there  are  relatively  large  increases 
in  combustion  efficiencies  under  low  operational  configurations.  We 
acknowledge  that  alterations  in  characteristic  length  are  key 

influences  on  the  degree  of  perfection  of  evaporation  or  vaporization 

* 

of  propellants.  From  results,  one  infers  that,  as  far  as  L  =0.5m 
variable  thrust  engines  are  concerned,  following  along  with  reductions 
in  thrust,  the  evaporation  or  vaporization  efficiency  is 
correspondingly  diminished.  This  is  in  line  with  the  results  of  the 
previous  theoretical  analysis.  ^ 

Fig. 8  gives  rj;  *,exp  values  for  L  =0.8m  and  L  =1.0m.  From 

the  comparison  in  the  Fig.,  it  is  possible  to  see  that  the  effects  on  'rj 

* 

* ,exp  are  not  obvious  when  L  rises  from  0.8m  to  1.0m. 
c  * 

Comparing  the  obvious  effects  on  tj  c*,exp  with  L  increased  from 

0.5m  to  1.0m,  and,  at  the  same  time,  considering  the  direct 

* 

responsivness  of  rj  *vap  to  L  ,  we  recognize  that,  speaking  in 

c  * 

terms  of  the  actual  engine  in  question,  when  L  =1.0m,  '  v  c*,vap 
already  tends  toward  a  saturated  value.  The  significance  of  this  is 
that  propellant  evaporation  or  vaporization,  as  far  as  the  various 
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operational  configurations  are  concerned,  is  already  basically 

complete.  It  is  no  longer  a  key  factor  influencing  efficiency  7 

* 

*,exp.  Because  of  this,  L  =1.0m's  7  *exp,  m  reality, 
c  .  c 

reflects  the  magnitudes  of  other  energy  losses  besides  mixing 
efficiency  and  evaporation  or  vaporization  losses. 


Fig. 8  A  Comparison  of  Efficiencies  With  Different  Combustion  Chamber 
Lengths  (II)  Solid  Symbol  L  0.8m  Hollow  Symbol  L  =1.0m  (1) 

Test 


Fig. 9  gives  actual  characteristic  speed  efficiencies  7,  *,exp 
★  C 
with  L  =  1.0m  and  the  cold  or  cool  flow  test  results  (  .7. 

*,mix)cool  for  the  same  jet  injection  device.  As  was  discussed 

c  * 

before,  when  L  =1.0m,  7  *,vap  ■»  1.0.  Because  of  this,  (  7 
*  *- 

*,exp)L  =1.0  =  7  *, mix, exp.  Comparing  (.7  *,mix)cool  and 
C  *  c  c 

(  ij  c*,exp)l  <  1.0m,  it  is  possible  to  know  that  the  difference  in 
values  between  them  follows  drops  in  operational  configuration  and 
increases.  Under  low  operational  configurations,  hot  test  results  (.17 
*  ,exp)  1  =  1 . 0m,  compared  to  cold  or  cool  test  results  7.  *mix,cool, 
will  be  much  lower.  The  explanation  for  this  is  that  cold  flow  test 
results  are  still  not  adequate  to  reflect  the  actual  mixing  processes 
of  the  engine  in  question.  At  the  same  time,  this  reminds  people 
that,  under  low  operational  configurations,  there  is  the  possibility 
of  the  existence  of  other  control  factors  related  to  energy  transfers. 
From  observations  during  engine  hot  tests,  following  along  with 
changes  in  operational  configurations,  there  are  clear  differences 
with  the  state  of  expelled  combustion  gases  in  jet  tubes.  Under  high 
operational  configurations,  the  combustion  gases  expelled  are  red  and 
bright.  However,  under  low  operational  configurations,  the  flame 
brightness  is  extremely  low.  Using  thermocouples  to  measure,  stagnent 
temperatures  were  approximately  1400°K.  This  will  be  much  lower 
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than  the  chemical  thermodynamic  equilibrium  temperature  for 
calculation,  2700°K,  as  compared  with  corresponding  conditions.  We 
know  that  there  is  a  close  relationship  between  the  speeds  at  which 
chemical  reactions  are  carried  out  and  temperatures  and  pressures. 
When  combustion  chamber  pressures  are  very  low,  the  non-equilibrium 
factors  associated  with  the  chemical  reactions  then  have  the 
possibility  of  becoming  one  of  the  control  factors  associated  with 
energy  transfers. 


Fig. 9  A  Comparison  of  Mixture  Efficiencies  in  Cool  and  Hot  Tests  (1) 
Cold  or  Cool  Test  Results  (2)  Hot  Test  Results 


V.  CONCLUSIONS 


1.  As  far  as  variable  thrust  liquid  rocket  motors  are  concerned, 
following  along  with  alterations  in  operational  configuration,  the 

combustion  efficiencies  will  show  the  occurence  of  changes.  With 

* 

regard  to  the  actual  L  =0.5m  engine,  when  the  thrust  is  adjusted 
down  in  a  5:1  alteration,  theory  and  experimentation  prove  that  the 
degree  of  imperfection  or  incompleteness  of  evaporation  (vaporization) 

and  mixing  follow  reductions  in  thrust  and  increase.  Tests  also 

* 

verify  that  appropriate  increases  in  the  length  L  are  capable  of 
causing  combustion  efficiencies  under  low  operational  configurations 
to  exhibit  relatively  wide  increases.  For  example,  an  engine  with 

it 

thrust  that  changes  from  125kg  to  25kg,  will  have  L  rise  from  0.5m 
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to  0.8m,  causing  combustion  efficiency  to  follow  the  range  of  the 
change  in  thrust  and  be  reduced  from  95%-60%  to  95%-75%.  Moreover 
the  increase  in  engine  weight  will  not  exceed  5%. 


r 


2.  Research  clearly  shows  that  the  key  factors  in  making 
propellant  mixing  processes  imperfect  or  incomplete  are  increases  in 
jet  injector  flow  speeds  during  injector  flow  adjustment  processes  and 
reductions  in  the  thickness  of  liquid  membranes  or  fluid  films. 

Making  fluid  films  or  liquid  membranes  unstable,  they  disintegrate  to 
become  strip  shaped  bundles  of  liquid,  causing  deviations  in  ram 
mixing  results.  Because  of  this,  in  designs,  it  is  necessary  to  make 
delicately  balanced  selections  for  the  jet  injector  degrees  of  opening 
h  and  h  .  within  the  flow  adjustment  range.  An  ideal  variaole 
thrust  liquid  rocket  motor  or  engine  will  execute  simultaneous  control 
of  propellant  component  ratios  and  jet  injection  speeds. 

3.  When  nitric  acid  27/asymmetrical  dimethyl  hydrazine 
propellants  are  used  in  variable  thrust  liquid  rocket  engines,  at  very 
low  pressures,  the  speeds  of  chemical  reactions  drop,  which  is  worthy 
of  serious  attention.  It  has  the  possibility  cf  becoming  one  among 
the  control  factors  for  energy  transfer. 
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